COMPUTING POLYNOMIALS OF THE RAMANUJAN t n CLASS 

INVARIANTS. 
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Abstract. We compute the minimal polynomials of the Ramanujan values 
t„, where n = 11 mod 24, using Shimura reciprocity law. These polynomials 
can be used for defining the Hilbert class field of the imaginary quadratic field 
Q(y , —n), and have much smaller coefficients than the Hilbert polynomials. 



Introduction 

Ramanujan on his third notebook, pages 392 and 393 in the pagination of |l( 
vol. 2] defined the values 



(1) 

where 



f(<ln /3 )f(ql) 



q n = cxp(-7rv / n). 



The function / is equal to: 

oo 

(2) f(-q):=l[(l- q n ) = q 



1/24 



n=l 



where q = exp(27rir), r g H and t](t) denotes the Dedekind eta- function. Without 
any further explanation on how he found them, Ramanujan gave the following table 
of polynomials p n (t) based on t n for five values of n: 



n 


Pn(t) 


11 


t - 1 


35 


t 2 + t - 1 


59 


t 3 + 2t - 1 


83 


t 3 + 2t 2 + 2t - 1 


107 


t 3 - 2t 2 + At - 1 



In Bruce C. Berndt and Heng Huat Chan proved that these polynomials indeed 
have roots the Ramnaujan values t n . Unfortunately, their method could not be 
applied for higher values of n and they asked for an efficient way of computing the 
polynomials p n for every n. Moreover, the authors proved that if the class number 
of K n := Q(V — n) is odd and n € N is squarefree so that n = 11 mod 24 then t n 
is a real unit generating the Hilbert class field. 

It is known that the Hilbert class field can also be constructed by considering 
the irreducible polynomial of the algebraic integer j(6) where 9 = — 1/2 + iy/n/2. 
The minimal polynomial of j{6) is called the Hilbert polynomial. It is interesting 
to point out that the coefficients of the polynomials p n have remarkably smaller 
size compared to the coefficients of the Hilbert polynomials. Therefore, finding 
an efficient and simple method for their construction is highly desirable and has 
a direct impact to applications where the explicit construction of class fields is 
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needed. Problems such as primality testing/proving |lj], the generation of elliptic 
curve parameters |7j and the representability of primes by quadratic forms Q can be 
considerably improved if the polynomials p n could be constructed with an efficient 
and easily implemented way. 

An explicit construction of the Hilbert class field has been done by N. Yui and D. 
Zagier in |l2] ] using the Weber functions. Yui and Zagier use a clever construction 
of a function on quadratic forms ax 2 + bxy + cy 2 that does not depend on the equiv- 
alence class of quadratic forms. The construction of a similar function in the case 
of p n polynomials seems very complicated and it is clear that a different approach 
must be followed. Our construction came from the enforcement of Shimura reci- 
procity law on the values t n . Shimura reciprocity law has been proven to be a very 
powerful tool for attacking similar problems g], ||, Q and can provide methods 
for systematically determining the instances when a given function yields a class 
invariant and for computing the minimum polynomial of a class invariant. 

The contribution of this paper is twofold. Firstly, we prove that the values t n 
constitute class invariants for all values of n = 11 mod 24. Expanding the theorem 
in Q , we show that t n generates the Hilbert class field not only in the case that the 
class number of K n is odd but also when it is even. Secondly we provide an efficient 
method for constructing the irreducible polynomials p n from the Ramanujan values 
t n and thus answer to the demand placed in Q for a direct and easily applicable 
construction method. Moreover, we have implemented our method in gp-pari J^] 
and we present all polynomials p n for all integers 107 < n < 1000, where n = 11 
mod 24 (table 1). 

The rest of the paper is organized as follows. In the first section we fix the 
notation and we give the Ramanujan t n values in terms of the Dcdekind eta func- 
tion. Next, we define six modular functions of level 72 i?, Ri, R2, R3, i?4, R§ and we 
compute the action of the generators of the group iSZ^Z) on them. In section 2 we 
prove that t n is indeed a class invariant, for all values n = 11 mod 24, and in final 
section we employ Shimura reciprocity law in order to compute the conjugates of 
t n under the action of the class group and compute the minimal polynomial of t n . 



1. Notation 

Let SLiiJl) be the group of matrices with integer entries and of determinant 
one. It is known [O, cor. 1.6] that the group SL^I) is generated by the matrices 

S= (J andT= (J I 

Every matrix ( a ^ j of induces an action on the upper half plane 



:= {z € C : Im(z) > 0}, 



by sending z >— » M+2 



cz-\-d ' 

Let 77 denote the Dedekind function: 



(3) t?(t) = exp(27rir/24) ]J (1 q n ), where r G H and q = exp(27rir). 

71=1 

Then 77-function is transformed by S and T as follows |ll], prop. 8.3] 

(4) 77(r + 1) = e 27rl/24 ?7(r) and ??(--) = V^?7(r). 

T 

From equation (Q) we can compute the action of every element g of 6X2(2) on the 
r;-function, since g can be written as a word in S*, T. 
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We denote by H n the Hubert field of K n := Q(-y/— n), i.e., the maximal Abelian 
unramified extension of K n . The extension H„/K n is Galois with Galois group 
equal to the class group of fractional ideals modulo principal fractional ideals. For 
imaginary quadratic fields the class group can be represented as the space of binary 
quadratic forms ax 2 + bxy + cy 2 modulo an equivalence relation ^, th. 5.30]. We 
will denote by [a, b, c] the quadratic form ax 2 + bxy + cy 2 and we will call two 
quadratic forms [a,,6j,Cj] for i = 1,2 equivalent if the corresponding roots t.; 6 H 
are in the same orbit of SL?($) acting on H. Using the identification of equivalence 
classes of quadratic forms with the ideal class group we can define the structure of 
an abelian group on the set of equivalence classes of quadratic forms. 

Let £q '■— (I7 1, (d = —n = 1 mod 4, n 6 N) be the zero element in this 
group. This element corresponds to the root 

1 , .v« 

Tl °=—2 +l — ■ 



Set 
Then 



q n = exp(-7n/n) = - exp(27riT£ ). 

/(?n) = /( - exp(27rir<> )) = exp(27riT £o )~ 1/24 f7(TA)), 
/(<£) = exp(27rzT, )- 3 / 24 7 ? (3r, ), 

Taking equation (Q) and all the above equations into consideration we arrive easily 
at the following: 

Lemma 1.1. The Ramanujan value t n is given by 

t n = V3i?2(r^ ), 

where 



r?(3T)?7(|r+ |) 
ry 2 (r) 



Now let N e N and r(JV) be the group 

r(iV) := | 7 e SL 2 (Z), 7 = ^ J M mod TV 

The field of modular functions of level TV consists of the meromorphic functions g of 
the upper half plane H that are invariant under the group r(iV), i.e. g{^r) — g(r) 
for every t G H and 7 E T(N). Every modular function is periodic with period N 
and thus it admits a Fourier expansion of the form 

00 

g(q) = X] avq ' ' ■> 

u=—i 

where q = exp(27rir/-/V). We will limit ourselves to modular functions where all 
coefficients of the Fourier expansions are elements of the field Q(Cjv)- The Galois 
group Gal(Q(0v)/Q) is isomorphic to the group (^) by defining a d (( N ) = 
for every (d,N) = 1. 

The action of the group Gal(Q(£jv)/Q) can be extended to the field of modular 
functions of level N with coefficients in Q(Gv), as follows: 



(5) g{q) ad = J2 °d{a v )q\ a d € Gal(Q(( N 



Moreover, the action of an element A £ GL2CZ) on modular functions g(q) can 
be expressed as 

(6) g A = [gBy^W 

where A = B ■ (l det(A)) and B £ S M^)- 

Lemma 1.2. The following functions are modular functions of level 72. 

R(r) - 

Mr) = 

Mr) = 
Mr) = 
R 4 (r) = 



V 2 (r) 




77(3r)r?(T/3 -} 


'1/3) 


rf(r) 




v(3t) V (t/3 -} 


'2/3) 


V 2 (r) 




V (r/3)v(r/3 - 


f 2/3) 


V 2 (r) 




77(r/3)77(r/3 - 


M/3) 



?y 2 (r) 

?7(t/3 + 2/3)^/3 + 1/3) 



t] 2 (t) 

Moreover, the element ad ■ C72 >— * C72 f or (d, n) — 1 acts on them as follows: 

(7) a d (R) = R, 

, . _ / C72 _1 ^i «/d=l mod 3 
<r*Wi) - I C d-2 R2 lfd ^ 2 mod3 



MR 



■2) 



( 2 2~ 2 R 2 ifd = l mod 3 
C72 _1 i?i i/d = 2 mod 3 

fr (n\-I ^2 _1 ^3 tfdsl mod 3 

^3j-| ^-2^ lfd ^ 2 mQd3 

, / C72 _2i ?4 i/d = 1 mod 3 

M ^ " 1 C?2 _1 i«3 tfd = 2 mod 3 

Proof. The fact that the above equations are indeed modular of level 72 is a direct 
computation using the transformations of the 77-functions under the generators T, S 
of 6X2(2) given in d). The action of given in (Q) is computed by considering 
the Fourier expansions of the 77-factors of the functions Ri. For instance let as 
compute the action of the element ad on R 2 . We begin by computing its action on 
77(r/3 + 2/3): 

(8) 77(7-/3 + 2/3) = cxp — ( T /3 + 2/3) £ a n exp( — -r + — -) = 



2^ = 

OO 



cxp (|£(t/3)) C 7 2 2 E Cr °« exp(^r), 



^=0 

*24 



where £3 = C72 is a primitive third root of unity. The desired formulas of equations 
(0) follow from the definition of the action of ad on (72 and arguing as above. For 
example (notice that from (p|) the Fourier expansion for rj(r) has rational coefficients 
so it is invariant under the action of a d ), 

°& (R2) = cr d jf-, = -YrT a dVn(r 3 + 2/3) . 

V V 2 (r) J V 2 {r) 



By (ft) we have that 



a d ( V (r/3 + 2/3)) = exp (^(r/3)) $ exp(^r). 



If d = 1 mod 3 then Q vd = Cs", thus 



2d-2 
72 



i?2. 



If c? = 2 mod 3 then £■ 



Q, and 



<7 d (77(r/3 + 2/3))=ry(T/3 + l/3))C 



2(2-1 



CTd(i?2 



/-2d-l p 
C, 72 rti. 



□ 



Later we will use some computer algebra programs in order to prove that t n is 
indeed a class invariant and find the minimal polynomials of t n . For this reason 
it is convenient to have the actions of the elements S, T, ad in matrix form. Using 
(0) we give the following matrix action of the elements 5*, T, ad on the functions 



R, Ri, 



(9) 



and 



/ R(r+1) \ 
Ri(t + 1) 
R 2 (t+1) 
R 3 (t + 1) 
Ra{t + 1) 
V fls(r + l) / 



.4 



/ R(r) \ 
Ri(r) 
R 2 (t) 
R 3 (r) 
R±{t) 
\ Mr) J 

( <J d R \ 

<?dRl 

OdR 2 
OdR-i 
OdRi 
V <JdR$ J 



= A 



( R(^) \ 
Ri(^) 
Rzi^r) 

R^) 
V M=±) J 

( R \ 

Ri 

i?2 

R 3 
Ri 

V ^ J 



= A f 



( R(r) \ 
Ri(t 
R 2 (t 
R 3 (t 
R a (t 
V Rs(t) J 



where 



(10) 



A T := 



/o 


C72 











°\ 








r 3 

S72 











r 6 

S72 





























1 

W2 




















1 


V° 








1 

<T2 








n 














°\ 











1 








C72 ( — C72 "^£72 ) 














C72 





-C 7 3 2 U + C7 B 2 





z-33 , a9 
S>72 T" C.72 




















^72 ^^72 











^72 


\o 














V 



An, — 



/I 

















\ 





/-d-1 
S>72 




















/■2d-d 
t >72 






















r 2d-2 
^72 






















/■d-1 
S>72 







\o 














<572 


7 



if d = 1 mod 3, 



A n , — 



/I 

















\ 








/■d-2 
S>72 














t >72 




























t 2d-\ 
S72 
















/■d-2 
<572 










Vo 














/-3d- 
S72 


7 



if d = 2 mod 3, 



2. SHIMURA RECIPROCITY LAW AND t„ CLASS INVARIANTS 

Let K n = Q{\/—n) be an imaginary quadratic number field, O the ring of 
integers of K n and = \ — i^-. It is known that j(9) is an algebraic integer that 
generates the Hilbert class field of K n , and moreover that the conjugates j(0) under 
the action of the class group are given by 



where T[ a ,b,c] is the unique root of ax 2 



= j{ T [a,b,c]), 

bx + c with positive imaginary part, i.e. 



„6,cl — 



2a 



There is an efficient algorithm for computing a set of non-equivalent quadratic 
forms and the minimal polynomial /d{x) of j(9) can easily be computed from the 
floating point approximation of the values j(t[ 0) (, )C ]): 

fo(x)= Yl ( x_ i{T[a,b,c])). 

[a,b,c]eCl(d) 

For instance, the polynomial for the quadratic extension of discriminant —107 is: 

/_ 107 (V) = a; 3 + 129783279616-10 3 a; 2 -6764523159552-10 6 x+337618789203968-10 9 . 

The disadvantage of using the above polynomials for the construction of the Hilbert 
class field is the very large size of their coefficients compared to the coefficients 
of the polynomials p n . Notice that the polynomial p n for n — 107 is equal to 
x 3 — 2x 2 + 4x — 1. In the literature there are alternative explicit constructions 
of the Hilbert class fields based on the Weber functions |[2]], or by other modular 
functions || , 0] , [pj . In [|| the authors proved that also the values t n can generate 
the Hilbert class field providing the following theorem. 

Theorem 2.1. If n = 11 mod 24 and the class group of K n is odd then t n gener- 
ates the Hilbert field. 

Proof. |, th. 4.1]. □ 
The Shimura reciprocity law can be applied in order to compute the minimal 



polynomial of t n and give an alternative proof of theorem 2.1 by removing the odd 
class number requirement. 

In order to prove that t n is a class invariant also for the case that the class group 
of K n is even, we will use the following construction: 

Theorem 2.2. Let O = Z[6>] be the ring of algebraic integers of the imaginary 
quadratic field K , and assume that x 2 + Bx + C is the minimal polynomial of 
9. Let N > 1 be a natural number, x%,...,x r be generators of the abelian group 
(O/NO)* andcti+Pid € O be a representative of the class of the generator Xi. For 
each representative we consider the matrix: 

i - Bfii -C3i 

Pi OLi 

If f is a modular function of level N and if for all matrices Ai it holds that 
(11) f{6) = f Ai {6), and Q(j) C 

6 



Ai 



then f(9) is a class invariant. 
Proof. |, Cor. 4] 



□ 



We know by lemma LI that the Ramanujan invariant t n can be constructed 
by evaluating the modular function y/3R,2 of level 72 at ti . Thus, we begin by 
constructing the generators of (0/720)* as Theorem 1.1 dictates. Since n = 11 
mod 24 we can take 

•44 

as a generator of the ring of algebraic integers of K = Q[y/ —n) . The minimal 
polynomial of 9 is x 2 — x + and thus J5 and C in Theorem 2.2 are equal to —1 
and ^j^- respectively. The form of the minimal polynomial implies that the prime 
p = 2 stays inert in the extension K/Q while the prime p = 3 splits. 

In order to prove that t n := ^/3R2(0) is indeed a class invariant we have to prove 
that 

(12) (\/SR2) A% = V3R2, for all matrices Ai and for n = 11 mod 24. 

We observe that the structure of the group (ffgy) depends only on the value of n 
mod 72, and there are exactly three equivalence classes n mod 72 so that n = 11 
mod 24, namely n — 11, 35, 59 mod 72. Thus, we have reduced ( |l2] ) to be checked 
for finite number of n. 

Using Chinese remainder theorem we can express the group (wj) as a direct 
product 

no) \§o) V80, 

We will study the structure of the above two summands separately. 
We compute that 

OV z z 
%0~) = 6Z X 6Z' 

A selection of generators for this group is cumbersome to do by hand. We have used 
brute force method, (we have checked all ellements one by one if there are invertible 
and then we have computed their orders) using magma |@j algebra system in order 
to compute that for C = G {3, 9, 15} a set of generators is given by 79 + 4, 5. 
Moreover the group 

o \* ^ z z z 

SO) = 12Z X 2Z X 2Z' 
and we have computed again using magma the following selection of generators 



n 


n+i 

4 


Generators 


11 


3 


6>,7,46' + 7 


35 


9 


59 + 6, 7, 46* + 7 


59 


15 


(9,7,461 + 7 



From the above generators and from the Chinese remainder theorem we can con- 
struct generators for the group (0/720)* and map them to the matrices Ai defined 
in theorem 2J2, We have totally 5 generators for the group (7^) : the first two 
are generators of the group (0/90)* and the last three are generators of the group 



(0/80)*. In order to compute the term / 1 in ([11]) we have to consider any lift of 
Ai in GL2(Z) and write it as a product Wi(S, T)diag(l, det(Aj)), where Wi(S,T) is 
a word in S,T. 

The following lemma gives us the decomposition of a matrix in 5X2(Z/p r Z) as a 
word in the generators of the group SX2(Z/p r Z). Therefore, we must consider the 
matrices Ai modulo 8 or 9 and define A^g € GL2(Z/8Z) and A^g £ GL2(Z/9Z), 
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so that Ai = Ai t s mod 8 and Ai = Ai t g mod 9. Notice that A^s = Id mod 8 for 
i = 1,2, i.e. the first two generators and Ai t g = Id mod 9 for i = 3,4,5, i.e. the 
last three generators. 

Lemma 2.3. Let p r be a prime power and let ^ € 5X2(Z/p r Z) so that 

either a or c is invertible modulo p r . Let S p r = f ^ Jj, T p r = ^ J J 6e too 

generators of the group SX 2 (Z/p r Z). Sei y = (1 + <x)c mod p r i/ (c,p) = 1, 
otherwise set z = (1 + c)a mod p r . T/ien 

a 6\ f T%rS p rT- r S p rTp- h mod p r (c, p) = 1 

c 1 S p J-^f/Vi; 4 r d modf/ tf(a,p) = l. 

Proof. (|, lemma 6]. □ 

The generators S p >-, T p r modulo p r = 8, 9 are then lifted to elements Ss , Sg , T% , Tg S 
SL 2 (Z) so that 

Ss = ( ^ mod 8, and Ss = Id mod 9 

Sg = ( ^ mod 9, and Sg = Id mod 8 

Tg = ^ ^\ mod 8, and Ts = Id mod 9 

Tg = ^ mod 9, and Tg = Id mod 8. 

Using Chinese remainder theorem we compute that 

S 8 = T^ST^ST^ST- 162 , T 8 = T 9 , 

- 1 Sg = T-^T-^ST-^T 1096 and T 9 = T~ 8 

We observe that the elements Ss,Tg commute with Sg,Tg modulo 72. 
The matrices Ai tP r , p r = 8 or 9 can be decomposed as products 

where the matrices Bi iP r have determinant 1 mod p r and can be expressed, using 
lemma 2.3, as words w p r (S, T) in the generators S, T. The matrices of the form Ai 
act on the field of modular functions of level 72 with coefficients in Q(C72) as the 
product of wg(S,T) ■ wg(S,T) ■ diag(l,dj), where di is the determinant of A i: i.e. 
the unique integer so that di = det^l^g) mod 9 and di = det(Ai^) mod 8. 

For instance the generator 70+4 of (0/90)* for = 3 corresponds to the matrix 

A := ^ \7 4/ mod 9 (and to the identity matrix modulo 8). This 

is a matrix of determinant 2 mod 9 and it is decomposed as 

'2 6\ _ (1 3\ fl N 
.7 4) ~ 1 7 2) U 2, 

/2 3\ . r— 1 

where B := I _, ^ ) i s a matrix of determinant 1 mod 9. Using we find that 

| 2 3\ m3n Ti7n t=i3 

— I 7 2 J ~~ J 9*->9Jg<39-f 9 - 

A lift of the elements Sg,Xg in SZ^Z) is given by equation (|l3|). This means that 
we replace each Sg of the above formula to Sg = t^ 1 ST~ 65 ST^ 1 ST 1096 and each 



Tg to T~ 9 . This gives us the desired lift of B to an element in SL^i^E). Using this 
lift, and the transformation matrices As, At given in ( |ic| ) we compute that the 
action of A on the modular functions Ri is given in terms of the following matrix: 

/ + \ 



E := 



c#-& 



^72 S72 





3 ' 3 



-Ct 2 o 

-2$l+e 7 2 



V o C 7 2+C 7 b 2 o ooo 

Let V be the Q( £72) -vector space of modular functions generated by the elements 
R, Rx, i?2, R3, R4, R5- The vector space V can be identified to the vector space 
Q(C72) 6 j m terms of the map 

V -> Q(C 72 ) 6 

a R + aiRi H h a 5 R 5 h-> (oo, Oi, ■ ■ ■ , a 5 ), 

a; e if om(V, Q(Ct2)). The space V* = Q(C 7 2) 6 is the dual space of V and we have 
to see how the action of the elements T,S,crd act on V*. The elements At, As 
defined in ( [Toj ) act on <Q(C 7 2) 6 in terms of the transpose matrices A T ,A S , while 
the action of A ad on R, Ri, . . . , R$ given on ( p"o| ) acts on Q(C 7 2) 6 in terms of the 
contragredient action, i.e. by considering the transpose of the matrix A a _ d . By 

Chinese remainder theorem we compute that the element i\ j? ] acts on Q(£ 7 2) 



v° 2 y 

as the automorphism ergs : (72 1— > C 7 |- Indeed, 65 is an integer 65 = 2 mod 9 and 
65 = 1 mod 8. 

Since d = 65 = 2 mod 3 we compute that the vector (0, 0, 1, 0, 0, 0)' correspond- 
ing to the element R 2 is mapped to A a _ 65 E a -^ (0, 0, 1, 0, 0, 0)* = (0,0,-1,0,0,0), 
where by _E CT - 65 we denote the matrix where all elements are acted on by cr_d. 
Notice that ^/i = C72 ~ C 7 2- Indeed, the value can be expressed as a differ- 
ence of two primitive 3-roots of unity £3, (f since i — C72 an d C3 = C 7 |- Moreover, 
(7-65 (\/3) = Cra 6 65 — Ct2 3 ° 65 — — V3 and thus -\Z3-R2 is left invariant. Following the 
same procedure it can be proven that \f3R2 stays invariant for all matrices A%. 

Theorem 2.4. The Ramanujan value t n is a class invariant for n = 11 mod 24. 



Proof. The condition Q(j) C Q(/) of theorem ||H is known j| proof of th. 4.1]. 



By machine computation 1 , it turns out that ( |l2| ) holds for all matrices Ai and thus 
t n is a class invariant for all n = 11 mod 24. □ 



3. Computing the polynomials p„ 

In this section we provide a method for the construction of the minimal poly- 
nomial of t n - Following the article of A. Gee ]5|, eq. 17] we give the following 
definition: 

Definition 3.1. Let N € N and [a, 6,c] be a representative of the equivalence class 
of an element in the class group. Let p be a prime number and p r be the maximum 
power of p that divides N . Assume that the discriminant D — b 2 — 4ac = 1 mod 4. 



The magma program used for this computation is available in the following web location 
http: //eloris . samos . aegean. gr /papers .html 
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We define the matrix 



a V 
1 

-6-1 



ifp\a 

— -c\ 
2 i/p I a tmdpfc 

if p \ a and p \ c. 



2 "2 
1 -1 



Chinese remainder theorem implies that 

GL 2 (Z/NZ) = Yl GL 2 (Z/p r Z). 

p\N 

We define ^4[ a ,b,c] as the unique element in GL 2 (Z/NZ) that it is mapped to A^ abc ^ p r 
modulo p r . This matrix A[ a b c ] can be written uniquely as a product 

(14) A [aM = B [aAc] ( I ^ ) , 

where d[ a ,f>,c] = det -Ar ,&,c] and Br a j )C ] * s a matrix with determinant 1. 

Shimura reciprocity law gives us |^|, lemma 20] the action of [a, b, c] on y/3R 2 (9) 
for 9 = 1/2- iVn/2: 

(\/3-R2(0)) ' a ' _6 ' C ' = (^ 6d [ a - 6 - c ] — £ 30d[o,I, ' cl ).R 2 ( Q [a,fe.c] T [a.fe.c] + ffia,i>,c] 

where | a [ a ' b ' c l A a >M J = A , . an d T \ a b c] is the (complex) root of az 2 + bz + c 

\7[o,6,c] °[a,b,c]J 1 ' ' ' 1 ' ' ' 

with positive imaginary part. 

If we try to implement this method in order to compute the polynomial for t n 
we face a problem. Even though we can compute a floating point approximation of 
the conjugate R 2 (A^ abc ^), it is n °t possible to use this approximation in order to 
compute the action of b , on it. There is however a simple approach that we 
can follow and solve this problem. We can express the matrix A[ abc \ as a product 
of a matrix Br a b c -i as in (^) and then compute the expansion of Bu bc \ as a word 
of the matrices S, T. 

We begin our computation by computing a full set of representatives of equiva- 
lence classes [a, b, c]. Since 72 = 2 3 • 3 2 we have to compute matrices 

A la , b<c] , p r G GL 2 (Z/p r Z) for p r = 8,9. 

Then we compute the determinant d[ a .&. c ],p r of the matrix A< a>biC i p r and we find a 
decomposition 

'\ 



A[ a ,b,c],p r — B[ a ^, c ],p r I n j 

The matrices •B[ 0l & )C ],p r are elements of SL 2 (Z/p r Z) and can be written as words of 
S p r,TpT by using lemma 2.3. 

So, if -Br a) j, )C i 8 = w(Ts,Ss) and Bt ab AQ — w'(Tg,Sg) are the decompositions of 
B[a,b,c],p r as words of S p r,T p r we take the lift 

B [a,b,c\ ^w(T 8 ,S 8 )w / (T 9 ,S 9 ) £ SL 2 (Z), 

and the corresponding action on the functions R, R\ 1 . . . , R$ is computed by using 

®- 

The determinants d[ a ,&, c ],8 € Z/8Z and d[ a ,f),c].9 € Z/9Z can be lifted in an 
element rfr a ,6,c] °f Z/72Z (so that it reduces to du^^g mod 8 and d\ a ,&, c ],9 mod 9 
respectively) by using again the Chinese remainder theorem. 
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The desired polynomial p n can then be computed: 

[a,b,c] ^ ' 

We have used the gp-pari 2 program in order to perform this computation. The 
resulting polynomials p n for 107 < n < 1000 are given in table 1. 
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Table 1. Polynomials p n for 107 < n < 1000. 



n 


Pn(t) 


107 


x 3 - 2x 2 + Ax - 1 




131 


x 5 + x A — x 3 — 3x 2 + 5x — 1 




155 


x 4 + 2x 3 + 5a; 2 + Ax - 1 




179 


x 5 - 2x 4 + 5a; 3 - x 2 + 6x - 1 




203 


x 4 - 3x 3 + 7x - 1 




227 


x 5 - 5x 4 + 9x 3 - 9x 2 + 9x - 1 




251 


x 7 + 5a; 6 + 6a; 5 - 2x 4 - Ax 3 + 2x 2 + 9x - 1 




275 


x 4 - x 3 + 6a; 2 - \\x + 1 




299 


x s + x 7 - x e - 12a; 5 + 16a; 4 - I2x 3 + Wx 2 - lSx + 1 




323 


x 4 - x 3 + Ax 2 + I3x - 1 




347 


x 5 + 7x 4 + 21x 3 + 27x 2 + I3x - 1 




371 


x 8 + 9x 6 - 10x 5 + Ux 4 + 8x 3 - 23x 2 + lSx - 1 




395 


x 8 - x 7 + Ux 6 + 16x 5 + 28x 4 + 2Ax 3 + 27x 2 + 17x - 1 




419 


x 9 - 6x 8 + 12x 7 - 7x 6 + I2x 5 ~ 8x 4 + 31x 3 + I0x 2 + 20x - 1 




443 


x 5 - 4a; 4 - 3x 3 + 17x 2 + 22x - 1 




467 


x 7 + 6x 6 + 7x 5 - 3a; 4 + 3x 3 - 23x 2 + 26x - 1 




491 


x 9 + x 8 + 16x 7 + 2x 6 + 37a; 5 - 31a; 4 + 44a; 3 - 40a; 2 + 29a; - 1 




515 


x 6 + 8x 5 + 32x 4 + 60a; 3 + 68a; 2 + 28a; - 1 




539 


x 8 - 6x 7 + 28a; 6 - 56a; 5 + 77a; 4 - 56a; 3 + 28a; 2 - 34a; + 1 




563 


x 9 + Ax 8 + 6x 7 - 11a; 6 + 44a; 5 - 76a; 4 + 91a; 3 - 64a; 2 + 38a; - 1 




587 


x 7 + x e + 16a; 5 - 12a; 4 + 20a; 3 + 24a; 2 + 39a; - 1 




611 


x w - 8x 9 + 35a; 8 - 62a; 7 - a; 6 + 116a; 5 - 65a; 4 - 100a; 3 + 125a; 2 - 


- 46a; + 1 


635 


a; 10 - 11a; 9 + 50a; 8 - 121a; 7 + 201a; 6 - 192.x 5 + 87a; 4 + 51a; 3 - 98a; 2 + 49.x - 1 


659 


x 11 - 7x 10 + 7x 9 + 27a; 8 + 19.t 7 - 43a; 6 - 5x 5 + 91.t 4 + 157a; 3 + 97.x 2 + 49a; - 1 


683 


a; 5 + 6a; 4 - 5x 3 - 41a; 2 + 56a; - 1 




707 


a; 6 + 4a; 5 + 30a; 4 + 72.x 3 + 108a; 2 + 58a; - 1 




731 


a; 12 + 7X 11 + 25a; 10 + 12a; 9 + 41a; 8 + 9a; 7 + 
+92a; 6 + 73a; 5 - 133a; 4 + 216a; 3 - 153a; 2 + 67a; - 1 




755 


a; 12 - 2X 11 + 18a; 10 + 50a; 9 + 82a; 8 + 182a; 7 + 360a; 6 + 522a; 5 + 
+598a; 4 + 486a; 3 + 262a; 2 + 66a; - 1 




779 


.a; 10 + 8x 9 + 24a; 8 - 8a; 7 - 11a; 6 + 26a; 5 + 81a; 4 + 220a; 3 + 98a; 2 + 


74a; - 1 


803 


a; 10 + 3x 9 + 26a; 8 + 11a; 7 - 65a; 6 + 16a; 5 + 7a; 4 - 83a; 3 + 150a; 2 - 


83a; + 1 


827 


a; 7 - 7.t 6 + 38a; 5 - 54a; 4 + 112a; 3 - 146a; 2 + 89a; - 1 




851 


a; 10 - 7x 9 - a; 8 + 86a; 7 + 69a; 6 - 201a; 5 - 219a; 4 + 94a; 3 + 103a; 2 - 


- 95a; + 1 


875 


a; 10 - 10a; 9 + 25a; 8 + 10a; 7 + 15.x 6 + 94a; 5 - 35a; 4 - 120a; 3 + 85a; 2 


+ 100a; - 1 


899 


a; 14 + 16a; 13 + 97a; 12 + 308a; 11 + 666a; 10 + 1086a; 9 + 






+ 1490x 8 + 1766a; 7 + 1800a; 6 + 1556a; 5 + 998a; 4 + 698a; 3 + 229a; 2 


+ 106a; - 1 


923 


a; 10 - a; 9 + 30a; 8 - 81a; 7 - 29a; 6 + 56a; 5 + 211a; 4 - 27a; 3 - 110a; 2 - 


- 115a; + 1 


947 


a; 5 + 5.x 4 + 7a; 3 - 103.x 2 + 125a; - 1 




971 


a; 15 - x 14 + 21a; 13 + 133a; 12 + 264a; 11 + 310a; 10 + 216a; 9 + 






+62a; 8 - 100a; 7 - 300a; 6 + 152.x 5 + 338a; 4 + 79a; 3 - 285a; 2 + 135a; 


- 1 


995 


.a; 8 + 12a; 7 + 59a; 6 + 78a; 5 + 12a; 4 + 66a; 3 + 289a; 2 + 140a; - 1 
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